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We present the theory of the fc-core pruning process (progressive removal of nodes with degree 
less than k) in uncorrelated random networks. We derive exact equations describing this process 
and the evolution of the network structure, and solve them numerically and, in the critical regime 
of the process, analytically. We show that the pruning process exhibits three different behaviors 
depending on whether the mean degree (q) of the initial network is above, equal to, or below the 
threshold {q)c corresponding to the emergence of the giant fc-core. We find that above the threshold 
the network relaxes exponentially to the fc-core. The system manifests the phenomenon known as 
“critical slowing down”, as the relaxation time diverges when (q) tends to (q)c. At the threshold, the 
dynamics become critical characterized by a power-law relaxation (oc 1/t^). Below the threshold, 
a long-lasting transient process (a “plateau” stage) occurs. This transient process ends with a 
collapse in which the entire network disappears completely. The duration of the process diverges 
when (q) —>■ {q)c. We show that the critical dynamics of the pruning are determined by branching 
processes of spreading damage. Clusters of nodes of degree exactly fc are the evolving substrate for 
these branching processes. Our theory completely describes this branching cascade of damage in 
uncorrelated networks by providing the time dependent distribution function of branching. These 
theoretical results are supported by our simulations of the fc-core pruning in Erdos-Renyi graphs. 

PACS numbers: 89.75.Fb, 64.60.aq, 05.70.Fh, 64.60.ah 


I. INTRODUCTION 

Pruning algorithms for networks provide an effective 
way to extract subgraphs distinguished by their struc¬ 
tural properties, connectivity, robustness against failures 
and damaging, and other features M- In general prun¬ 
ing processes, parts of a network are progressively re¬ 
moved from it according to some rule. If the rule is sim¬ 
ply random removal of nodes, we obtain ordinary per¬ 
colation but in general we are interested in more 

complex pruning rules. The parts removed may be nodes 
[IIjSOj clusters m , finite connected c omp onents in in¬ 
terdependent and multiplex networks [6I. II2 h 14| . etc. De¬ 
spite the wide variety of pruning processes, many of them 
demonstrate similar behaviors, such as discontinuous hy¬ 
brid phase transitions. The fc-core pruning as the sim¬ 
plest pruning process of this kind, stands as a paradigm 
for all such pruning processes, so its theory should help 
to understand the behavior of these pruning algorithms 
in general. The fc-core is the network subgraph in which 
all nodes have degree at least fc [ 1 ]. Since fc-cores repre¬ 
sent the densest parts of networks, they play an impor¬ 
tant role in understanding the structure and dynamics 
of complex network systems Q . The standard algorithm 
for finding the fc-core of a network employs the following 
pruning process: at each step remove all nodes of degree 
less than fc. This removal decreases the degrees of re¬ 
maining nodes, some of which will become smaller than 
fc. So, the pruning is repeated until either the fc-core 
remains or the network disappears [a. 
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Previous investigations have mainly focused on the fi¬ 
nal result of the fc-core pruning process, namely the fc- 
core. These were the studies which showed that fc-core 
percolation is a hybrid phase transition, combining dis¬ 
continuity and a critical singularity, in contrast to ordi¬ 
nary percolation (continuous phase transition) 0, [l^ . 
However, associating the number of steps in the pruning 
process with time t reveals a process exhibiting complex 
dynamics above, below, and at the fc-core percolation 
threshold. Understanding the fc-core pruning process and 
accompanying structural changes can shed light on such 
physical phenomena as the jamming transition, the rigid¬ 
ity percolation, and glassy dynamics [l^. Furthermore, 
the fc-core pruning process is one of the simplest exam¬ 
ples of dynamic processes associated with hybrid phase 
transitions, sharing, for example, some common prop¬ 
erties with cascade failures in interdependent networks 
that have recently received significant attention in the 
literature IT^. [l3l. [TtM^ . 

In this paper we develop the detailed theory of the fc- 
core pruning process in uncorrelated, sparse random net¬ 
works, describing the temporal evolution of the networks 
structure, the spreading of damage over the network, and 
critical phenomena in this process. We show that near 
the threshold value of the mean degree, {q)c, correspond¬ 
ing to the emergence of the giant fc-core, this cascade of 
removals of nodes is a branching process with the mean 
branching coefficient close to 1. Our theory describes 
this process completely providing the full time depen¬ 
dent distribution of branching from the beginning until 
the end of the pruning. We indicate that the clusters of 
nodes of degree fc (so-called “corona clusters”), evolving 
due to the pruning, provide the substrate for the branch- 
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ing processes. Near the threshold we find three different 
behaviors depending on whether the mean degree (q) of 
the initial network is above, equal to, or below the thresh¬ 
old. First, we demonstrate that above the threshold, the 
network relaxes exponentially to the steady k-core. The 
relaxation time diverges when {q) tends to {q)c, mani¬ 
festing a phenomenon known as “critical slowing down”. 
Second, at the critical point, (q) = {q)c, the dynamics 
is critical, characterizied by a power-law relaxation with 
l/t^ dependence. Third, below the threshold, a long- 
lasting transient process (a “plateau” stage) occurs. This 
transient process ends with a collapse in which the entire 
network disappears. We find that the duration of the 
process diverges when (q) approaches {q)c- Our theory 
is supported by numerical calculations for Erdos-Renyi 
graphs and by simulations of the pruning process in these 
random graphs. 

In Sec. m we derive the exact equations describing the 
evolution of the network structure during the pruning 
process enabling us to obtain the time dependent degree 
distribution P(q,t) and the branching probability distri¬ 
bution P(n, t) at all times. Close to the critical point, the 
probability that different branches of the process cross 
each other is negiligibly small. We show that in this 
region, our equations take a simple form for analytical 
treatment. Section Hill explores the three regimes of the 
pruning process below, at, and above the threshold. Sec¬ 
tion m describes the statistics of the branching process. 
A relationship with dynamical systems close to a saddle- 
point bifurcation and details of calculations are given in 
the appendices. 


II. EVOLUTION EQUATIONS 


To study the k-core pruning process, let us consider as 
a representative case an infinite uncorrelated sparse ran¬ 
dom network, which is completely defined by its degree 
distribution P{q). In this case, we can write exact equa¬ 
tions for the evolution of the degree distribution. Let 
P{q,t) be the proportion of vertices having degree q at 
time t, with the initial condition P{q, 0) = P{q). At each 
time t = 1,2,3,..., all vertices with degree q less than 
k are pruned by having all edges connected to them re¬ 
moved from the network. The probability P{0, t) thus 
tracks the number of vertices pruned so far. 


A. Exact evolution equations 


The removal of edges from pruned vertices means that 
some non-pruned vertices will also lose edges, changing 
the degree distribution of the remaining network. Let 
Tt be the probability that, upon following a randomly 
chosen edge within the network existing at time t, we 


arrive at a vertex with degree less than k: 

U = (1) 

q<k 

Such an edge will be removed in the subsequent step. 
Here {q)t is the mean degree of the surviving network at 
time t, 


{q)t = ^qP{q,t). ( 2 ) 

q 

The probability that a vertex of degree q' > k a.t 
time t has q surviving edges at time t -I- 1 is then 
. A vertex of degree q' < k dX time t 
will of course have degree zero at time t 1. Summing 
over all q', the degree distribution then evolves as follows: 

P{q,t + l)= P{q' (3) 

q' 

for g > 0 while the fraction of pruned nodes evolves ac¬ 
cording to 


P(0,t + 1)= 


( 4 ) 


q' <k 


where the sum includes q' = Q. The uncorrelated na¬ 
ture of the network ensures that Eqs. CD"® completely 
define the evolution of the network at all times. Note 
that another approach for the pruning process which, 
however, does not consider the evolution of the network 


structure, was proposed in |22|. 

To understand the spreading of damage through the 
network as the pruning process evolves, we introduce the 
probability St- This is the probability that, following an 
edge at time t, we reach a vertex that has degree at least k 
at time t, but will have no more than k—1 other surviving 
edges at time < -I- 1 (not counting the edge through which 
we reached the vertex). This means that if the edge we 
are following is removed at time t, the vertex that it leads 
to will be removed at time t -I- 1. To calculate St, we sum 
over probabilities that all but I of the q—1 outgoing edges 
of a vertex of degree q (i.e. q — 1 — I edges) are lost at 
time t (each one with probability rt) with I equal to at 
most A: — 1. A second summation is then performed over 
all possible degrees q > k: 


St 


q>k 1=0 ^ ' 


The probability P{n,t) that a vertex removed at time t 
has n neighbors that will be removed at time t -|-1 is then 


V{n,t) = 


Y!rJnP{<l.t){l)sni-StY- 

Y!rJrP{<l,t) 


( 6 ) 
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FIG. 1. (Color online) A snapshot of the branching process of propagation of node pruning in a small part of the network of 
10® nodes during the plateau stage ((g) < (g)c) of the fc-core pruning process for fc = 3. The node labelled 0 is pruned, causing 
the corona nodes (i.e., nodes with degree 3) labelled 1 to lose edges. These two nodes are pruned in the next step, and so on, 
with further corona nodes removed in subsequent steps according to the numbered order. Red and white circles represent the 
nearest neighboring nodes of degree 4 and greater than 4, respectively, that survive because their degrees exceed 3. The red 
nodes after this pruning become of degree 3. They augment other corona clusters, which may then be pruned at a later time. 


This function describes the branching of spreading dam¬ 
age. The mean branching is 


fc-i 

5t = ^ n'P{n,t) 

n—0 


EqJ 

EqZlPiq,t)' 


(7) 


B. Non-crossing approximation 


Unfortunately, it is difficult to study analytically 
Eqs. (Pi-dll). In this subsection we develop an approx¬ 
imate approach providing the asymptotic description of 
the pruning process at large times near the critical point. 

When the probability r* is very small, the pruning 
can then be considered as a branching process. The 
probability that a vertex loses two neighbors in a sin¬ 
gle step is negligible, in other words, the probability 
that two or more branching trees meet at a vertex is 
negligible. The process then evolves with independent 
branching trees spreading simultaneously over the net¬ 
work. An example of such non-crossing branchings ob¬ 
served in simulations is shown in Fig. P If crossings 
are negligible, then the fraction of vertices of degree 
q < k — 1 is also negligible and only vertices of degree 
q > k — 1 must be taken into account. This is the main 
assumption of the “non-crossing” approximation. This 
approximation is supported by our numerical solution of 
Eqs. dU^dU simulations which show that the prob¬ 
ability of crossings between branches are negligible and 
P{k — 1, t) P{k — 2,t) » ... P(l, t) already after a 
short initial period (see the next sections). Applying the 
non-crossing approximation to Eq. (0, we find that St 
becomes simply the probability that, following an edge 
at time t, we encounter a vertex of degree k. 


St 


kP{k, t) 
{Q)t 


( 8 ) 


Furthermore, the probability r^, Eq. dm, and the mean 
branching bt, Eq. d3, take the simple forms, 


n 


bt 


{k — l)P{k — 1, t) 

{<l)t 

(k — l)kP(k, t) 

{Q)t 


(9) 

( 10 ) 


So rt is simply the probability that, following an edge 
at time t, we encounter a vertex of degree k — 1. The 
evolution equation (0 is also simplified. The following 
set of equations determines the evolution of the degree 
distribution during the k-core pruning process: 


P{q, t+l)=P{q, t)-rtqP{q, t)+rt{q+l)P{q+l, t), (11) 
P{k — l,t + l)=rtkP{k,t), (12) 

P(0,t-t-l) = P(0,t)-fP(/c-l,t), (13) 


{q)t = (k- l)P(k + qP{q, t), (14) 

q>k 

where q > k. The negative term in Eq. m corresponds 
to the reduction in P{q,t) due to vertices of degree q 
losing with the probability qrt a single edge, while the 
positive term (last term) corresponds to an increase in 
P(g, t) due to vertices of degree q -f 1 losing an edge with 
the probability (g -f l)rt and so ending up with degree q. 
Using Eq. m, we rewrite Eq. dni) as follows, 

P{k-l,t + l) = btPik-l,t). (15) 

Equations m and m show that the removal of a ver¬ 
tex of degree fc—1 at time t triggers in the next step the 
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removal of all corona vertices attached to it since they 
will lose one edge and will have degree k—\. On av¬ 
erage, the number of these corona vertices is the mean 
branching ht- In uncorrelated networks, Eqs. ([TT1)-(fn|) 
describe the non-crossing branching processes of spread¬ 
ing damage (see Appendices iBl and ICll. They show that 
vertices of degree k (“corona” vertices) are crucial for 
spreading damage. In the case (q) > {q)c at large times, 
t ^ 1, crossings are negligible and these equations are 
asymptotically exact. Equations (HHi-dll are not valid 
when there are numerous crossings between branching 
processes. Such crossings are abundant both at the ini¬ 
tial stage of the pruning process and at the end of the 
“plateau” stage when the network collapses. In this case, 
the exact Eqs. 0-0 must be used. Branching processes 
are discussed in detail below in Sec. El 


III. THREE REGIMES OF THE PRUNING 
PROGESS 

As a representative example of the pruning process, we 
solved Eqs. 0-0 numerically for Erdos-Renyi networks 
(Poisson degree distributions) using the initial mean de¬ 
gree (g) as a control parameter. We solved the equations 
for fc = 3 and fc = 5. The 3-core appears with a hy¬ 
brid transition at {q)c ~ 3.35091887, while for the 5-core, 
(g)c ~ 6.7992755. We also performed simulations of the 
pruning process in the networks. We found that for any 
fc > 3, the dynamics of the pruning process can be di¬ 
vided into three different regimes: (g) < (g)c, (g) > (g)c, 
and (g) = (g)c- 


A. Pruning process below (g)c 

Below (g)c, the pruning process ends in a finite time 
(number of steps) with the complete destruction of the 
infinite network. Rapid pruning of vertices at early times 
soon slows down and the system enters a “plateau” stage 
in which the rate of removal of vertices is very slow. Fi¬ 
nally, this transient process ends with a collapse in which 
the entire network disappears, as can be seen in Fig. 
which displays the temporal dependence of the network 
size S. The duration of the entire process, from begin¬ 
ning until final collapse, diverges as the inverse square 
root of the distance from the critical point, 

T = Abeiow/\/(g)c - (g), (16) 

as shown in Fig. |3l The time T is mainly determined by 
the duration of the “plateau” stage. Note that the in¬ 
verse square-root scaling law is a general feature of non¬ 
linear dynamic systems that are close to a saddle-node 
bifurcation [H, [^. In such systems, the long-lasting 
transient process is caused by a bottleneck region (the 
ghost) that exists in phase space when the system is 
close to a saddle-node bifurcation or the limiting point 


of metastable states in the case of the first order phase 
transitions (see a simple model in Appendix lAl). The na¬ 
ture of the bottleneck effect in the fc-core pruning process 
is discussed in Sec. E 



FIG. 2. (Color online) Size S of the Erdos-Renyi network vs 
time t during the pruning process for fc = 3 in two cases: (i) 
below the threshold {q)c, the system passes through a long 
“plateau” stage before a final collapse. Shown are numerical 
calculations for mean degree (g) = 3.3509 (blue solid line) and 
simulations (triangles) for a network of 10® vertices showing 
similar total time, (ii) Above (g)c, the system relaxes to a 
finite size, numerical solution for (q) = 3.35092 (red solid 
line) and simulations (circles). 


B. Pruning process above (g)c 

Above (g)c a finite fraction of the network remains in¬ 
definitely and the network relaxes to the steady fc-core 
only in the infinite time limit (see Fig.[2|). In this regime, 
according to the numerical solution of Eqs. o-® and 
simulations, the relaxation to the steady state is exponen¬ 
tial. Instead of measuring the total time, we characterize 
the time scale of the pruning process by measuring the 
relaxation time r, where: 

P{k — l,t) (x . (17) 

The relaxation time r diverges as the inverse square root 
of the distance from the critical point, as seen in Fig. [31 

— Aabove/\/ { q } c (g) • (1^) 

We examine the origin of this scaling in more detail in the 
next Sec. lIVBi using the non-crossing approximation. 

The divergence of r manifests the phenomenon known 
as critical slowing down. Furthermore, comparing the 
amplitudes Abeiow and Aabove of the square-root singu¬ 
larities below and above the transition, we find their ratio 
to be Abeiow/Aabove = 9.133/1.452 = 6.29 for fc = 3 and 
8.44/1.34 = 6.28 for fc = 5, in agreement with the ra¬ 
tio 27r expected for general transitions of this kind, see 
Appendix [^ 
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FIG. 3. (Color online) Characteristic times associated with 
the k-core pruning process for fc = 3 on an Erdos-Renyi net¬ 
work. Circles show the duration T of the entire pruning pro¬ 
cess below (q)c. Squares show the relaxation time constant r 
above {q)c- Fitted square root scaling for T and r are shown 
by black solid lines. The critical point {q)c is marked by a 
vertical dotted line. Completion times for a simulated net¬ 
work of 10® vertices are also shown below {q)c (triangles). In 
the inset, the inverse squares of T and r (also scaled by 27r) 
are shown, demonstrating the inverse square dependence on 
the distance from the critical point. 


In Fig. 0] we show the evolution of the degree distri¬ 
bution just above and just below {q)c. Near the critical 
mean degree the initial evolution of the degree distri¬ 
bution P{q,t) both above and below the critical point is 
similar, namely, there is a sharp initial decrease of P{q, t) 
for nonzero q. Below {q)c, however, the network finally 
collapses completely, while above the critical point, the 
fc-core survives forever. The theoretical results agree well 
with simulation. 


C. Critical pruning process 

Solving Eqs. 0-0 numerically for Erdds-Renyi net¬ 
works, we find that exactly at the critical point, {q)c, the 
relaxation is much slower, with P{k — 1, t) decaying as a 
power law, 

P(/c-l,t)cx^. (19) 

as can be seen in Fig. [5] For fc = 3 we measured the 
exponent a = —1.993 at {q) = 3.35091887, suggesting 
that the exponent is —2. Note that in a simple model of 
a particle moving in a one-dimensional potential in Ap¬ 
pendix]^ the corresponding critical exponent is —1 [see 
Eq. mi We explain the power-law behavior, Eq. (HU, 
in Appendix iBlbv solving Eqs. (IB1I) - (IB4I) within the non¬ 
crossing approximation. This approach gives the exact 
value a = 2. 


IV. BRANCHING PROCESSES OF PRUNING 

In this section, to understand the nature of the critical 
dynamics of the pruning process, we study the spread¬ 
ing of damage through the network and the structural 
changes during this process. The probability P(n,t) of 
the branching process is given by Eq. m that takes a 
simple form within the non-crossing approximation, 

P(n,t)= (20) 

Since the parameter st given by Eq. ([5]) is the probabil¬ 
ity to encounter a vertex of degree k (corona vertex), the 
probability Eq. (l20t is precisely the probability that fol¬ 
lowing an edge we arrive at a corona vertex which has n 
corona vertices at the ends of emanating edges [l^. It 
is important to note that as the network evolves accord¬ 
ing Eqs. 0-0 during the k-core pruning process, so 
do the corona clusters and hence their size distribution. 
Since the probability s* and, therefore, the probability 
V(n, t) of the branching process depend on time, the size 
distribution of branches of removed vertices is therefore 
related but not equal to the instantaneous size distribu¬ 
tion of corona clusters (see the following subsection). 

A. Branching processes at (q) < {q)c 

The numerical solution of Eqs. 0-0 and simulation 
show that, in the case (q) < {q)c during the plateau stage, 
the pruning process develops as a branching process, as 
described in Sec. Ill Bl The branching process of removals 
evolves in agreement with Eqs. (O-dlSl). Examples of 
typical pruning trees are illustrated in Fig. [5] for k = 3. 
The figure shows that crossings between the branching 
trees are only abundant at the beginning of the pruning 
process and are rare in the plateau stage. The crossings 
also are abundant at the end of the plateau stage, signal¬ 
ing a collapse in which the entire network disappears. 

The full branching distribution, given by Eq. m is 
shown in Fig. [T] It is similar both above and below the 
transition, and barely changes throughout most of the 
pruning process. In the figure we also show the branch¬ 
ing distribution observed in simulations. The agreement 
with theory is good, however there are noticeable finite 
size fluctuations, which are largest when the pruning pro¬ 
cess is slowest: this occurs in the middle of the plateau 
period. In contrast to this behavior, fluctuations in the 
case {q) > {q)c are enhanced with increasing time (see 
Fig. 7). Critical behavior of fluctuations is a common 
property of systems approaching the critical point of a 
continuous phase transition, or the limiting point of the 
metastable states of a first-order phase transition, how¬ 
ever, discussion of these phenomena is beyond the scope 
of the present paper. 

In the early part of the plateau, P{k — l,t) decreases, 
reaching a minimum in the middle of the plateau stage. 
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FIG. 4. (Color online) Time evolution of the network degree distribution during the fc-core pruning process, for an Erdos-Renyi 
network with k = 3. Each line shows P{q, t) for a different value of q, in order from top to bottom, g = 3,4, 5, 2, 1, as labelled, 
(left) Initial mean degree (q) = 3.3509, (right) initial mean degree (q) = 3.35092. Also shown are traces from simulation runs for 
an Erdos-Renyi network with N = 10® vertices and mean degree 3.3511 (left) and 3.35111 (right). Note that the critical point 
for a particular realisation is a stochastic quantity, so the mean degrees for matching theory and simulation are not necessarily 
equal. Theoretical curves were chosen to be near the critical point and to have a similar total time. 



FIG. 5. Decay of P{k — 1) on a log-log scale for several values 
of (g) close to {q)c- Exactly at the critical point, the decay 
follows a power law with exponent —2. (For clarity, the final 
collapse of P{k — 1) for curves below the critical point is not 
shown.) Curves are plotted for the Erdos-Renyi graphs with 
the mean degree values (labelled 1 to 6) 3.350905, 3.35091, 
3.350915, which are below q^, 3.35091887 (very close to gc), 
and 3.35092, 3.350925, which are above gc. Dashed line is a 
power-law decay with exponent —2. 


at t = tjn- From Eq. (IB2I) this corresponds to the point 
when bt reaches 1. A Taylor expansion around this point 
(see Appendix [C]) gives the temporal behavior of P{k — 
l,t) in the plateau stage 


The corresponding equation for bt is 

bt^l+(l- ^)Cb, ( 22 ) 

Our analysis of the plateau st age in Appendix ICl shows 
that Ch ^ IjT oc \/{q)c — (q) "C 1. This analytical re¬ 
sult agrees with our observation from the numerical solu¬ 
tions. The mean branching bt, Eq. 0, is slightly below 1 
in the beginning of the plateau stage. As time increases, 
bt tends to increase, as pruning of vertices decreases the 
mean degree of the network. The mean branching reaches 
1 at tm as required, then continues to increase, with an 
accelerating rate of pruning, until the network finally col¬ 
lapses rapidly, as seen on the left in Fig. HI We observe 
from numerical solution of the exact Eqs. 0-0, and 
from simulations, that the minimum occurs in the mid¬ 
dle of the plateau stage, i.e. tm = T/2, see Fig. H) 
The numerical solution of exact Eqs. 0-0 shows that 
the coefficient Cp is of order 1. Using this result and 
Eq. (IC5|) in Appendix [0 we find a relationship between 
P{k — l,tm) and the plateau duration T, 

P(fc-1,U) - ^ oc (g)c-(g). (23) 

The instantaneous size distribution 11(5', t) of finite 
corona clusters can be found directly from the degree 
distribution P{q,t) at every time t: 

U{S,t) = (24) 


P{k-l,t) = Pik-l,t^) 




, ( 21 ) 


where S*{t) —> cx) at the critical point of the emergence 
of a giant corona cluster. According to [1^, a giant con- 
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time 


FIG. 6. (Color online) Example of the pruning process from time t = 10 till 34 in a small part of the Erdos-Renyi network with 
10® nodes. Time progresses from left to right in the tree. Blue circles represent vertices removed at a given time step. Their 
removal results in the removal of vertices on the right, and so on. Crossings (dashed lines) between the branching processes 
are abundant at the beginning of the pruning process. They appear rarely after a short initial period. Typical trees of medium 
size are shown. Much longer and much shorter trees also occur. 
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FIG. 7. (Golor online) Evolution of the branching distribution Vijiyt) below (left) and above (right) the critical point for fc = 3. 
Solid and dashed curves are theoretical curves from Eq. (|6l) for n as labelled, for {q) = 3.3509 (left) and (g) = 3.35092 (right). 
Points are measured from simulation of an Erdos-Renyi network of 10® nodes, at (g) = 3.3511 (left) and 3.35111 (right), being 
just below and just above the critical mean degree for that network. 


nected cluster of corona vertices is present when 
{k-l)kP{k,t) 

bt = -TT- > 1. 

Wt 


(25) 


In the case of fc = 3, we have S*{t) = — l/ln[4s^(l — St)] 
where St = 3P{3,t)/{q)t according to Eq. (|S]) [l^. In 
Appendix Owe show that a giant corona cluster appears 
continuously at the same time tm when the fraction P{k— 
l,t) of fc — 1 nodes achieves a minimum. Such a giant 
corona cluster will be consumed by the pruning process, 
guaranteeing the collapse of the whole network in finite 
time. A similar continuous emergence of a giant subgraph 


prone to failure was recently observed in interdependent 
networks in Ref. [l^ . The left side of Eq. (USD is identical 
to Eq. (fTOl) . so the border of the region where a giant 
corona cluster appears is at the point where the mean 
branching of the pruning process equals 1. The region 
in the (g) — t plane where the giant corona cluster is 
present is marked in Eig.|9]as region III. Note that a giant 
corona cluster only appears below {q)c in the plateau 
stage. At (g) = (g)c, at any time t there are only finite 
corona clusters. When t —^ oo, the size distribution of 
corona clusters tends to the power law function Eq. (24), 
corresponding to the critical point of the emergence of a 
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FIG. 8. (Color online) Evolution of the mean branching ratio bt below (left) and above (right) the critical point for k = 3. 
Dashed curves are calculated using Eq. cni, points are from simulations. Parameters used and simulation realisations are the 
same as in Fig. [T] 



(5) 


FIG. 9. (Color online) Phase diagram for the fc-core pruning 
process in {q) — t plane. The vertical line represent the critical 
point {q) = {q)c (fc = 3 for this figure). In region II at (q) > 
{q)c, the pruning process reduces the network to the giant k- 
core as time approaches infinity. Only finite corona clusters 
are present in region II. A giant corona cluster is present in 
region III. The mean branching is 1 on the border between 
regions II and III. The mean branching is below 1 in region II 
and larger than 1 in region III. The network collapses at times 
on the upper boundary of region III so there is no network in 
region I. 


giant corona cluster. Above {q)c, there are only finite 
corona clusters at any time. 


B. Branching processes at (q) > {q)c 

Above the transition point, with increasing time the 
degree distribution P{q,t) tends to the steady distribu¬ 
tion Pk{q) with mean degree {q)k = J2q>k^^k{q) while 
P{k — 1, t) —> 0. In turn, the mean branching bt saturates 
at a constant value bk less than I (see the right side of 
Fig. IS]). If (g) is close to {q)c, I - k{k - l)Pk{k)/{q)k 


B\/{q)c — {q) where i? is a constant [ 2 ^, and using 
Eq. m we have ~ 1 — {q)c — (g). Substituting 
the constant bk for bt = k{k— l)Pk{k,t)/{q)t in Eq. (IB2I1 
in Appendix[Bl and solving, we find an exponential decay 
of P{k — l,t), Eq. (fTTl) . and a relationship between the 
relaxation time r and the branching coefficient bk, 

bk = l-T~^. (26) 

Therefore, 

= 1 - k{k - l)Pk{k)/{q)k ~ B^/ (g)^ - (g) (27) 

in agreement with the numerical solution Eq. (1181) . The 
pruning process only evolves within finite corona clusters, 
and the network survives at any time t (the region II in 
Fig. in]) approaching the steady A:-core as time approaches 
infinity. 

C. Critical branching process 

Exactly at the critical point, (g) = (q)c, the branch¬ 
ing bt comes arbitrarily close to 1, but only reaches that 
value in the infinite time limit. The leading term in 
1 — is a monotonically decreasing function of t. Solving 
Eqs. (IBll) - (jB2l) in Appendix iBl we find that the function 
P{k—1, t) has power-law behavior, Eq. (IT^ . with critical 
exponent a = 2. This behavior corresponds to the mean 
branching bt increasing as 

bt = l-2/t + 0{l/f). (28) 

This kind of time dependence of the mean branching 
is known to lead to the avalanche lifetime distribution 
L{T) oc found in various models (see, for ex- 

ample . 12 71 . l28[) and real systems (for example, in the 
brain [29j). This suggests that the power-law relaxation 
Eq. (IT^ and the avalanche lifetime distribution have the 
same origin. 
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Since the mean branching bt = (k — l)st tends to 1 
when t —)> oo, we have s* —>■ l/{k — 1). Equation (1^ 
gives the following exact result: 


P(n, oo) 


fk-l\ (fc- 2)'=-!-'^ 

\ n ) {k — 


(29) 


In the case of k = 3, we obtain P(0, oo) = P(2,oo) = 
1/4 and P(l,oo) = 1/2. These values agree with results 
obtained by our simulations and numerical solutions that 
are displayed in Fig. [T] 


V. DISCUSSIONS AND CONCLUSIONS 

In this paper we have developed the theory of the k- 
core pruning process in uncorrelated, sparse random net¬ 
works. Employing the numerical solution of the exact 
evolution equations, Eqs. 0-0, an asymptotic analy¬ 
sis, and simulations in Erdos-Renyi graphs, we revealed 
that this process demonstrates three different kinds of 
critical behavior depending on whether the mean degree 
{q) of the initial network is above, equal to, or below the 
critical point, {q)c, corresponding to the emergence of 
the giant fc-core. We found that above the critical point, 
at large times the network relaxes exponentially to the 
steady fc-core. At the critical point, (q) = {q)c, the dy¬ 
namics is critical and it is described by a power law with 
respect to time (oc 1/t^). Below the critical point, the 
pruning eliminates an infinite network in a finite time. 
The duration of the transient process diverges when {q) 
tends to {q)c from below. 

We found mechanisms for these critical phenomena. 
Studying the structure of paths along which damage is 
spreading in the network, we found that the damage 
spreading is a branching process. Our analysis showed 
that it is the evolving clusters of nodes of degree k 
(“corona clusters”) that provide the substrate for the 
branching process. Indeed, if a vertex of degree k loses 
an edge and is removed, then this removal triggers a re¬ 
moval of all corona vertices, one by one, which belong 
to the same corona cluster. Using analytical methods 
and simulation, we showed that the pruning can be con¬ 
sidered as a branching process that begins after a short 
initial period of rapid network change. During this pro¬ 
cess, independent branching trees develop with branch¬ 
ing ratio close to I. The temporal behavior of the mean 
branching plays a crucial role in the branching process 
and slowdown of the fc-core pruning dynamics at the crit¬ 
ical point and during the plateau stage. To understand 
the branching process it is important to note that corona 
clusters evolve in time. When damage propagates over 
the network, on one hand, it removes corona nodes, but 
on the other hand, it decreases the degrees of neighboring 
nodes, producing new corona vertices and thus increas¬ 
ing size of other corona clusters, which can be pruned 
at a later time. Due to this, the branching probability 
becomes time dependent. The mean branching is close 
to 1 during the whole plateau stage, below {q)c. At the 


beginning of the stage the mean branching is a little bit 
smaller than I, and slowly increases with time. It reaches 
the value I approximately at the middle of this stage and 
then continues to increase. At this point a giant corona 
cluster is formed providing a substrate for the complete 
collapse of the network at the end of the plateau stage. 
Exactly at the critical point, the mean branching comes 
arbitrarily close to I, but never reaches it at any finite 
time. This leads to a power law decay in the fraction 
of nodes of degree fc — 1. The branching trees of prun¬ 
ing become arbitrarily long, but a giant corona cluster is 
never formed, until t —>■ oo. Finally, we found that above 
the critical point {q)c, the mean branching saturates at 
a constant value less than I. In this case, mean size of 
branches is finite and relaxation to the steady fc-core fol¬ 
lows an exponential law. 

The fc-core pruning process in sparse, uncorrelated 
random complex networks is a representative model of 
dynamics in complex systems undergoing hybrid phase 
transitions. We have solved this model and have de¬ 
veloped the complete description of critical dynamical 
phenomena including the long-lasting transient process, 
critical relaxation, and critical slowing down. We sug¬ 
gest that our results could be useful for understanding 
similar collective phenomena that occur in other com¬ 
plex systems near discontinuous (hybrid and first-order) 
phase transitions. 
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Appendix A: Relaxation in ID system near border 
of metastability 


The behavior described in this paper for the fc-core 
pruning process is common for dynamical systems having 
a saddle point for some set of system parameters. Here 
we consider the simplest dynamical model of this sort, 
namely a particle moving in a one-dimensional potential 
F(x) in a viscous medium, demonstrating features similar 
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FIG. 10. Potential F{x) versus x in Eq. (lAlll in the following 
cases: a) the coefficient a > Uc = , F'{x) < 0 at any x. b) 

a = ac, F'{x = b) = 0 and there is a saddle at x = b. c) 
a < ac, F{x) has a local minimnm and a local maximum. 


Here r is the relaxation time, 


1 

2-\/5^ — a 


(AS) 


diverging at the critical point. 

The square root critical singularities of T and r, 
Eqs. and (IMI) . respectively, coincide with those of 
the k-coie pruning process, Eqs. (HU) and (HU. From 
these expressions we obtain the remarkably beautiful ra¬ 
tio of the critical amplitudes of T and r: 


to the k-coTe pruning process: 

dtx = -d^F{x), 

F{x) = —ax + bx^ — ia:^. (Al) 

o 

Here the coefficients a and b are positive, and the variable 
(particle’s coordinate) x{t) > 0. The initial condition 
is x{t = 0) = 0. There are three distinct regimes, see 
Fig. [TUI 

(a) a > Qc = b^ , normal phase, with dxF(x) < 0 at any 
X. At the end of the process x approaches infinity; 

(b) a = Oc, resulting in the saddle point Xg = bhi F(x). 
As t ^ oo, x(t) approaches b; 

(c) a < Qc = b^, which gives the local minimum 

(“metastable state”) at Xm = b — — a. As 

t —>■ oo, x{t) approaches Xm- 


The straightforward solution of Eq. dni shows that 
in the normal phase [regime (a)], i.e., when a > Oc, the 
variable x, starting from 0, approaches infinity in a finite 
time T{a). For a close to ac, the process greatly slows 
down when x{t) passes the value b, and we obtain the 
asymptotic expression 


a — 2bx + x'^ y/a - 6^ ' 

This time diverges at the critical point Oc = b^ [regime 
(b)], at which x relaxes to the saddle point value b slowly, 
in a power-law way. Asymptotically, we get 

x{t) -b^ (A3) 

Note that in the case F{x) = cx'^, which corresponds 
to the second order phase transition within the Landau 
theory, the equation dtX = —dxF{x) leads to the critical 
relaxation x oc . 

In regime (c), x(t) relaxes exponentially to the local 
minimum value Xm = b — y/b'^ — a. Asymptotically, 

x{t) — Xm OC (A4) 



T{a-ac) 

r(ac-a) 


(A6) 


coinciding with the corresponding ratio that we found for 
the fc-core pruning. 

Near the critical point, the time T to complete the 
process (to run away from x = 0 to infinity) in the nor¬ 
mal phase is strongly influenced by small variations of 
the parameters of the system. To quantify this effect, 
we introduce a time dependent perturbation h(t) of the 
coefficient a in the potential F{x), namely, a — h(t). Let 
h{t) be a constant h within the interval of width e around 
some moment i. We define the response of T to h as 


X{i) = 


lim 
h ,€—^0 


AT(/i, e, i) 
he 


(A7) 


where AT{h, e, i) = T{h, e, i) — T{h=0, e=0, i) is the vari¬ 
ation of the time T due to the perturbation h{t). This re¬ 
sponse takes an elegant asymptotic form as a approaches 
the critical point ac = b^, 



which has a Lorentz shape in terms of the moment of 
the perturbation, t, and diverges according to the Curie- 
Weiss law. Here T is the time, given by Eq. (IA2I) . to run 
away to infinity in the absence of perturbation, h = 0, 
and T /2 is the time at which the particle passes the point 
X = b. This divergence of the response x f^e critical 
point indicates the presence of strong fluctuations near 
the critical point, which we observe in the k-core pruning 
process (see Figs. Hand [5]). 


Appendix B: Critical relaxation in the non-crossing 
approximation 

Let us solve Eqs. (fni)^(|TT|) and find the critical behav¬ 
ior of P{k — 1, t) at the critical point {q) = {q)c- At t ^ 1, 
we consider P{q,t) as a function of continuous time t. In 
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this limit, Eqs. (inii-dii take a differential form, 

{q+l)P{q+l,t)-qP{q,t) (,B1) 


dP{k-l,t) 

dP{0,t) 


fc(fc-l) 


dt 


. {<l)t 

=P{k-l,t), 


P{k,t) - 1 


P{k-l,t), (B2) 


(B3) 


((?)t=(fc - l)P{k - 1, t) + ^ qP{q, t), (B4) 

q>k 

where q > k. In order to solve these equations, we use the 
fact that with increasing t the degree distribution P{q, t) 
ioi q > k tends to the steady degree distribution of the 
fc-core, Pk{q), i.e., Piq,t) = Pk{q) +SP{q,t). Moreover, 
6P{q, t) —>■ 0 and P{k—1, f) —>■ 0 in the infinite time limit. 
At the critical point, the distribution satisfies the condi¬ 
tion k{k - l)Pfc(fc) = {q)k where {q)k = Y.q>klPk{q) 
sm. We solve Eq. (IBII) in the first order in P[k — 1, t) 
and find 


5P{q,t)= 


{k-l) 


{q) 


k ^ 


qPk{q)-iq+l)Pk{q+l) JP{k-l,t)dt. 

(B5) 

Then, using Eq. (IB4|) . we find {q)t in the first order in 
P{k — Substituting these results into Eq. (jB2l) gives 
an equation, 

dP{k-l,t) 


dt 


/ OO 

P{k — 1, t)dt 

P^{k-l,t)+0{P^k-l,t)), (B6) 


{q)h 


where 


fc(fc—1)^(A:-|-1) 


Pk{k+l)- 


m 

Equation (ESI has a solution 


{k-l){k-2) 

{q)k 


P{k-l,t) = — +0{l/t^). 


(B7) 


(B8) 


Numerical solution of the exact Eqs. 0-0 confirms this 
result (see Eig.[5]). Using Eqs. (IB5|) and (IB8I) . we find the 
mean branching, 

{q)t t 

Note that this is the universal critical behavior of branch¬ 
ing processes [^ . 


Appendix C: Plateau stage in the non-crossing 
approximation 

If {q)<{q)c, with increasing time the fraction P(fc—1, t) 
of nodes of degree k — 1 achieves a minimum at a time 


tm in the middle of the plateau stage (see Fig. E). The 
time tjn is determined by the condition 


dPik-l,t) 


dt 


= 0 . 


(Cl) 


According to Eq. 
also holds. 


at t — t jYi the following equality 

k{k - l)P{k,tm) 


{q)t 


= 1 . 


(C2) 


It signals the percolation of corona clusters [see Eq. (EH)]. 
Thus, the minimum of P(fc — 1, t) occurs when the giant 
corona cluster appears. Near the minimum, we can use 
the Taylor expansion 

P(A:-I,t) =P(fc-l,t^)[l + i(^-l)"c'p + ...], (C3) 
where 

tl, d^P{k-l,t) 


P(fc-l,tm) 


dP 


(C4) 


Differentiating Eq. (IB2I) with respect to t, we find the 
second derivative and 


Cp - t^VuiP(,k 1, 


(C5) 


where 


'^m — 


k{k — I)^(fc 4- I) 
~Wr. 


P(k -\- 1, tm) ~ 


{k-l){k-2) 
{q)t 


(C6) 

We estimate P{k — l,tm) and tm using the numerical 
solution of exact Eqs. (HD-®. Our numerical results in 
Secs. IIII Al and HV Al show that the coefficient Cp is of the 
order of 1, and tm ~ P/2, i.e., the minimum takes place 
at the middle of the plateau stage. Equation (IC5|) gives 
a relationship between P(fc — l,tm) and the duration T 
of the entire pruning process, 

P(fc - 1, U) - ^ oc {q)c - (q). (C7) 

Note that in the neighborhood of the threshold {q)c, the 
plateau duration tends to the time T to complete the 
pruning process. Equation m and the Taylor expansion 
of the function P{k, t) give the temporal behavior of the 
mean branching bt near tm, 


bt = l + 






(C8) 


where 


Cb=tm ^^\^ ^'’ P{k-l,tm) 

\^)tm 

X \{k-lfP{k + l,tm)P2{k-l)-{q)tJ^. (C9) 

This equation shows that Cb is small since C{, ^ 1/P oc 
'J{q)c — (q) ^ 1- This result is also supported by our nu¬ 
merical solution and simulations for Erdos-Renyi graphs 
(see Sec. II V A]) . 
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